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HYERS ULAM STABILITY OF HIGHER-ORDER CAUCHY-EULER DYNAMIC 

EQUATIONS ON TIME SCALES 

DOUGLAS R. ANDERSON 

Abstract. We establish the stability of higher-order linear nonhomogeneous Cauchy-Euler dynamic equa- 
tions on time scales in the sense of Hyers and Ulam. That is, if an approximate solution of a higher-order 
Cauchy-Euler equation exists, then there exists an exact solution to that dynamic equation that is close to 
the approximate one. Some examples illustrate the applicability of the main results. 



< 

U; 1. INTRODUCTION 

In 1940, Ulam [26] posed the following problem concerning the stability of functional equations: give 
conditions in order for a linear mapping near an approximately linear mapping to exist. The problem 
for the case of approximately additive mappings was solved by Hyers [TT] who proved that the Cauchy 
equation is stable in Banach spaces, and the result of Hyers was generalized by Rassias [23J. Alsina and 
Ger [TJ were the first authors who investigated the Hyers-Ulam stability of a differential equation. 

Since then there has been a significant amount of interest in Hyers-Ulam stability, especially in relation 
to ordinary differential equations, for example see [3 El H21 [131 HH US H3 [201 E21 [25]. Also of interest 
are many of the articles in a special issue guest edited by Rassias [23], dealing with Ulam, Hyers-Ulam, 
\ and Hyers-Ulam- Rassias stability in various contexts. Also see Li and Shen [13 [18], Wang, Zhou, and 
Sun [27], and Popa et al [SUE]- Andras and Meszaros [3] recently used an operator approach to show the 

b ■ 

stability of linear dynamic equations on time scales with constant coefficients, as well as for certain integral 
equations. Anderson et al [21 Corollary 2.6] proved the following concerning second-order non-homogeneous 
Cauchy-Euler equations on time scales: 

Theorem 1.1 (Cauchy-Euler Equation). Let Ai, A2 € R (or A2 = Ai, the complex conjugate) be such that 

t + X k fj,(t) ^ 0, fc = l,2 
for all t £ [a, a{b)]f, where a G T satisfies a > 0. Then the Cauchy-Euler equation 

* AA (t) + 1 ~*Iq X2 * A (t) + ^| *(*) = f(t), t € [a, % (1.1) 
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has Hyers-Ulam stability on [a,b]j. To wit, if there exists y G C^[a, b]j that satisfies 

< e 



y (*) + v (*) + W) y{t) - f(t) 



for t G [a, 6]t, ^en i/iere exists a solution u G C^j 2 [a, 6]t o/ (|l.lj> given £>?/ 

u(t) = ex x (t,T 2 )y (t 2 ) + ex x (t,a(s))w(s)As, any t 2 <E [a,a 2 (b)]i 



where for any t\ G [a, ct(6)]t t^e function w is given by 



w{s) = ex 2 -i (s,ti) 
CT ( S ) 



y A (n) - ^y(ri)l + r e ^( S ,a(C))/(C)AC 



suc/i i/iai |y — u[ < .Ke on [a, a 2 (b)]j for some constant K > 0. 

The motivation for this work is to extend Theorem 1 1.1 1 to the general nth-order Cauchy-Euler dynamic 
equation; we will show the stability in the sense of Hyers and Ulam of the equation 

n 

J2a k M k y(t) = f(t), M y(t) := y(t), M k+1 y(t) := tp(t) (M k y) A (t), k = 0, 1, • • • , n - 1. 

fc=0 

This is essentially [5, (2.14)] if <^>(i) = t and /(t) = 0. Throughout this work we assume the reader has a 
working knowledge of time scales as can be found in Bohner and Peterson [HE], originally introduced by 
Hilger [10]. 



2. Hyers-Ulam stability for higher-order Cauchy-Euler dynamic equations 

In this section we establish the Hyers-Ulam stability of the higher-order non-homogeneous Cauchy-Euler 
dynamic equation on time scales of the form 



a k M k y(t) = f(t), M y(t) := y(t), M k+1 y(t) := <p(t) (M k y) A (t), k = 0, 1, 



, n — 1 



(2.1) 



k=0 



for given constants ajGl with a n = 1, and for functions <p, f G C r d[a, b]j, using the following definition. 

Definition 2.1 (Hyers-Ulam stability). Let (p,f G C r d[a, b]j and n G N. If whenever M k x G C^[a, 6]t 
satisfies 



^a k M k x(t) -f{t) 



k=0 



< e, t G [a, b]i 



there exists a solution it of (|2.ip with M/%u G CA[a, o]t for k = 0, 1, • • • , n — 1 such that |x — u| < on 
[a, a n (6)]T for some constant K > 0, then (|2.ip has Hyers-Ulam stability [a, b]j. 
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Remark 2.2. Before proving the Hyers-Ulam stability of (|2.1|) we will need the following lemma, which 
allows us to factor (|2.ip using the elementary symmetric polynomials [7\ in the n symbols pi, ■ ■ ■ ,p n given 
by 



»i = si(pi,--- ,p n ) = Y^ Pi 

i 

«2 = S 2 {pi,--- ,p n ) = ^PiPj 



i<j 



4 = S 3 (pi, ■ ■ ■ , p n ) = ^ pipjpk 



i<j<k 



si = S 4 (pi, ■ • • ,p n ) = ^ PiPjPkPl 

i<j<k<£ 



StiPl, ■ ■ ■ ,Pn) = PhPi2---Pk 

h<i2<—<h 



Sn = S n (pi, ■■■ ,p n ) = p 1 p 2 p3 ■■■ p n - 

In general, we let represent the ith elementary symmetric polynomial on j symbols. Then, given the 
in (I2.ip . introduce the characteristic values Xk G C via the elementary symmetric polynomial on the n 
symbols — Ai, • • • , — A n , where 

Oik = Sn-k = Sn-k( — \i, ■ ■ ■ , -A n ) = (- \) n ~ k \ h X i2 • • • X in _ k , a n = s = 1. (2.2) 

h<i2 <---<i„-k 

Lemma 2.3 (Factorization). Given y,tp G C r d[a, b]j and otk G M with a n = 1, Ze£ M^y G C~j[a, b]i, where 
Moy(t) := y(t) and Mk+iy(t) := y?(i) (M^y) A (t) /or fc = 0, 1, • • • , n — 1. XTien we Ziaue £/te factorization 

n n 

J2 a k M k y(t) = Y[(<pD- X k I) y(t), n G N, (2.3) 
fc=o fc=l 

where the differential operator D is defined via Dx = x A for x G C^[a, b]j, and I is the identity operator. 
Proof. We proceed by mathematical induction on n G N, utilizing the substitution defined in (|2.2p . For 

71 = 1, 

n 

£ a kM k y{t) = a M y{t) + ai M iy (t) = si(-Ai)l/(*) + 1 • ^(t)y A (t) = (pi) - AiJ) 

fc=0 
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and the result holds. Assume fj2.3j) holds for n > 1. Then we have a n+ \ = 1 and 

n+l n 

^a fc M fc y(t) = a y(i) + ^ afcMfcy W + M ™+ iy W 

k=l 

n 

+ E 'SU-WMrt*) + V(<) (M n y) A (t) 
fc=l 

n 

+ £ ( s n+i-k ~ ^y(t) + (M n y) (t) 



fc=0 fc=l 

ft 



Vi+1 



fc=l 

71 



k=l 



+ J2 C + i- k M k y(t) + <p(t)D (M n y) (t) 



k=l 

n 



-Xn+l s n-kM k y(t) + <p(t)D £ <+i-k M k-iV(t) + M n y (t) 

k=0 \k=l / 

n /n—1 \ 
-Xn+l s n-kMkV(t) + <p(t)D £ s n-k M ky(t) + M n y (t) 



k=0 \k=0 

n n 



-X n+1 s n n - k M k y(t) + <p(t)D £ aS_ fc Mjfcy(t) 

fc=0 fc=0 
n 

( V (t)D - \ n+l I)J2s^ k M k y(t) 

k=0 

n 

(<p(t)D - \ n+1 I)^2a k M k y(t) 

k=0 

n 

(<p(t)D - K+iI) ]J (<pD - A fe J) y(t) 



fc=i 

and the proof is complete. □ 

Theorem 2.4 (Hyers-Ulam Stability). Given y,<p,f € C r( j[a, b]j with \ip\ > A > for some constant A, 
and a k G M witt a n = 1, consider (|2.ip mtft M^y € C^j[a, o]t /or fe = 0, • • • , n — 1. Using the X k from the 
factorization in Lemma \2li\, assume 

<p(t) + \ h n(t) + 0, fc = l J 2,---,n (2.4) 

for all t € [a, o" n ~ 1 (6)]x- Then (|2.1I) /ias Hyers-Ulam stability on [a,b]j. 

Proof. Let e > be given, and suppose there is a function x, with M/^x € C^j[a, 6]t, that satisfies 



J2a k M k x(t)- f{t) 



k=0 



<e, te [a,b]j. 
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We will show there exists a solution u of (|2.ip with M^u G C^j[a, 6]t for k 
\x — u\ < Ke on [a, o" n (6)]T for some constant K > 0. 
To this end, set 

gi = (^x A — Ai x = (ipD — Ail) x 
52 = <pgi ~ A 2 ^i = (<pD - A 2 J) 5i 



0, 1, • • • , n — 1 such that 



5fc 



This implies by Lemma 12.31 that 



VSn-l ~ *n5n-l = - A n J) n _i. 



0»(t) - /(*) = ^a fe M fc x(t) - /(i), 

k=0 



so that 

M*)-/(t)| <e, t€[o,6] T . 
By the construction of g n we have Ipg^-i — A n g n _i — / | < e, that is 



A f 

9n-l 9n-l 

if if 



e e 
< — < — . 

- m - A 



By [21 Lemma 2.3] and (|2.4p there exists a solution u)i G C^j[a, b]j of 

ip(t) ip(t) 
t G [a, 6]t, where wi is given by 

ft ff s ) 

wi(t) = exn(t,Ti)g n -i(Ti) + eA 2L (t,cr(s))— T As, any n € [a, a(b)]r 

V J T1 V <P{8) 

and there exists an Li > such that 

\g n -i(t) < £ie/A, t € [a,<r(6)] T . 

Since #„_i = 9?5^_ 2 - K-i9n-2, we have that 

\f9n-2 -K-l9n-2(t) ~ Wl{t)\ < L\E j A, t G [a,cr(b)] T . 
Again we apply [21 Lemma 2.3] to see that there exists a solution u) 2 G C A [a, a(b)]j of 



(2.5) 
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t G [a, ct(6)]t, where W2 is given by 



w 2 (t) = e V i {t,T 2 )gn-2(r 2 ) + / e Vl (t, cr(s))— -— As, any r 2 G [a, cr 2 (6)] T , 
and there exists an L 2 > such that 

|s„- 2 (t) - w 2 (t)\ < LzL-te/A 2 , t G [a, a 2 (6)] T . 
Continuing in this manner, we see that for k = 1, 2, ■ ■ ■ , n — 1 there exists a solution -u^ G C^Ja, a k ~ 1 (b)]f 

of 

/(t)-^^)-^ = 0, or <p(t)vA(t)-\ n - k+1 w(t)-w k - 1 (t) = 0, 
(p(t) tp(t) 

t G [a,a h ~ l (b)]j, where w k is given by 

pt ( \ 

w k (t) = e x n _ k+l (t, r fc )ff n _ fc (r fc ) + / e A n _ fc+1 (t,o-(s)) fc ~\ As, any r k £ [a,a k (b)] T , (2.6) 

— 7 Tfc — <p{s) 

and there exists an L k > such that 

k 



\g n -k(t) - w k (t)\ <Y[L j£ /A k , t G [a,a k (b)] T . 

3=1 



In particular, for k = n — 1, 



n— 1 



| 5 i(t) -w„_i(t) | < IJ^M n_1 . * G K^" _1 (6)]t 



implies by the definition of g± that 

n-l 



<JjL je /A n , *€ [a,a n - 1 (&)]T. 
i=i 



Thus there exists a solution w n G C^[a, cr n 1 (6)]t of 

0, or tp(t)w A (t) - Aiio(t) - io n _i(t) = 0, 



U> (t) T-^-W[t) 



<p(t) w 

t G [a, (j ra ~ 1 (5)]ir, where w n is given by 

rt 



w n {t) = e\ 1 {t,T n )x(T n )+ I ex 1 (t,a{s)) As, any r n G [a, cj™(6)]t, (2.7) 
and there exists an L n > such that 

n 

|x(t) -«;„(*) | < Ke :=\{L ] e/A n , t G [a,a n (b)] T . (2.8) 

j'=i 
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By construction, 

(ipD - Ail) w n (t) = w n -i(t) 

2 

Y[ (pD - A fe J) w n (t) = (ipD - X 2 I) w n -i(t) = io n _ 2 (i) 

k=l 

n 

H(<pD-\ k I)w n (t) = (tpD - X n I) wi(t) = f(t) 

k=l 

on [a, a n ^ 1 (b)]j, so that u = w n is a solution of (|2.1j) . with u G C^j[a, cr n_1 (6)]T and |a:(i) — itf n (*)| < ^£ 
for t € [a, o" n (&)]ir by (j2.8j) . Moreover, using (|2.7j) and (|2.6p . we have an iterative formula for this solution 
« = w n in terms of the function x given at the beginning of the proof. □ 
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